GENERALIZED JACOBI OPERATORS IN KREIN SPACES 
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Abstract. A special class of generalized Jacobi operators which are self-adjoint in Krein 
spaces is presented. A description of the resolvent set of such operators in terms of solutions 
of the corresponding recurrence relations is given. In particular, special attention is paid to 
the periodic generalized Jacobi operators. Finally, the spectral properties of generalized Jacobi 
operators are applied to prove convergence results for Pade approximants. 
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Let /i be a positive Borel measure having infinite support supp a C [a, b] C R. To every 
such a measure « there corresponds a linear functional defined on the linear space V = span{A fc : 
k eZ + :=NU {0}} by the formula 

s k = &(\ k ) ■ = [ t k dfi(t), keZ+:=NU {0}. 

J a 

The functional 6 is positive definite on V, that is, det(sj +J )" J=0 > for all n G Z + . Besides, 
the measure « (or, equivalently, the functional &) generates the holomorphic function 



cn 

where R is large enough. By using the Euclidean algorithm, P.L. Tchebyshev [31] expanded 
the function /i into the following continued fraction 

, , 1 1 bl b\ 

(i.2) m - 2 — = — 



x % A - a - A - d - A - a 2 



A - a -J 



A — (Li — 

where dj are real numbers, bj are positive numbers, and sup{|aj| + bj} < oo. Such continued 



fractions are called J-fractions fH\. Note that the coefficients aj and bj are uniquely determined 
by the coefficients Sj of the Taylor series at infinity (see (11. ip ). 

It is well-known (see [5], [21], [28]), that the n-th convergent — Q n (A)/P n (A) of the con- 
tinued fraction fll.2p is characterized by the following property 

n q \ Qn(A) S Si S 2 „-i n ( 1 A v 

(L3) -pj\)= — -v ^ + {x^) (A ^ oo) " 

In other words, the rational function —Q n jP n is the n-th diagonal Pade approximant to /x, that 

is, 
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Further, by using standard argumentation (see [21]), one can see that the polynomials P n and 
Q n are solutions of the three-term recurrence relations 

(1.4) bj-iUj_i + cijUj + bjU j+ i = Xuj (j e N), 
with initial conditions 

(1.5) Po(A) = l, P 1 (A) = ^, g (A) = o, gx(A) = i 

On the other hand, to the recurrence relations (11.4j) (or, equivalently, to the continued 
fraction (jl.2p ) there corresponds a linear bounded self-adjoint operator in the space n. More 
precisely, that operator is generated by the following tridiagonal matrix 



H 



( a b \ 
b o x bx 



b\ a 2 



\ 



and its m-function m(A) = ((if — A) 1 e, e) £2 , where e = (1, 0, . . . , 0, . . . ) T G i? coincides 

[O.oo) I- 1 / 

with ju(A). So, one has an operator representation of // 

^(A) = ((H-A)- 1 e,e) <a , A e C \ R. 

[0,oo) 

The relations between orthogonal polynomials, Pade approximants, and Jacobi operators 
are well-known. These relations allow to use operator methods to the investigation of orthogonal 
polynomials and the Pade approximants. 

Note, that the above-mentioned results are also valid if the underlying functional (5 is 
not positive and has the property 

(1.6) det(s i+j )l j=0 ? 0, Si = 6(A') 

for all n e Z + (see PQ, [6], [7]). The functional & having the property (11.61) is called regular. 

In the present paper a similar relations for non-regular functionals are presented. In 
fact, the scheme proposed in PQ, [7] for investigation of the convergence of Pade approximants 
is generalized here to the non-regular case. The paper is organized as follows. In Section 2, 
starting from a (not necessarily regular) functional on V, three-term recurrence relations for 
associated polynomials are derived. In Section 3, a special class of generalized Jacobi matrices 
is presented. Weyl solutions of the three-term recurrence relations and Weyl functions are 
introduced in Section 4. In Section 5, following the scheme proposed in p], the characterization 
of resolvent sets of generalized Jacobi operators is obtained. Section 6 is concerned with the case 
of periodic generalized Jacobi matrices. In Section 7, convergence results for Pade approximants 
are proved. 

2. Associated polynomials 

In this section, starting from a (not necessarily regular) functional on V, three-term 
recurrence relations for associated polynomials are derived. 

Let us consider a holomorphic in a neighborhood of infinity function (p such that 



p»(A) := <p(A) = <p(\). 
So, <p has the Taylor expansion at infinity 



GENERALIZED JACOBI OPERATORS IN KREIN SPACES 



3 



where Sj G R. To every such a function one can associate a real linear functional on V defined 
by the formula 

6(A fc ) -.= —& X k ip(X)dX = s k e E, fceZ + :=NU {0}, 
2tt? J| A | =iJ 

for sufficiently large i?. Clearly, the functional & is not necessarily regular, that is det(sj +J -)" J - =0 = 
may vanish for some n G Z + (for instance, see [13]). In general case, the functional & gener- 
ates an indefinite inner product on V (see [32], [3]) 

[/, g] e := 6(/(A)<?»(A)) = -L / /(A)^(A)^(A)dA, f,g £ V, 

2m J\x\ =R 

which is degenerate if and only if is rational. In what follows we suppose that ip is not 
rational. As in the regular case, one can associate to & the following holomorphic function 

(2 ' 1} 62 v^Aj = 2^i N/? T3X = -A-p A^L"- (|A|>/?) ' 

Throughout this paper we suppose that the sequence s := {sj}°51 is normalized, i.e. 

the first nonvanishing moment has modulus 1. A number rij G N is called a normal index if 

det(sj + fc)^ =0 7^ 0. Since ip is not rational, there exists an infinite number of normal indices 

(see [TIJ Section 16.10.2]). Let n\ < n 2 < ■ ■ ■ < rij < . . . be a sequence of all normal indices. 

By the choice of n\ one has s m _i ^ 0. Let us set Eq = s m _i (|eo| = 1) and ipo := <p. The 

1 

principal part of the Laurent expansion for is a polynomial of degree ko := ri\ with the 

leading coefficient e . So, we have 

(2.2) = e po(X) + fe(A), 6 > 

<MA) 

where is a monic polynomial of degree ko and (pi is holomorphic in a neighborhood of infinity. 
Furthermore, the function <p 1 satisfies the relation <p\(X) = <Px(X)- Choose b > such that the 
sequence sW = {sj }^ defined by the following expansion 

(i) (i) „(i) 

(2.3) ^(A)= S ° Sl 2{n ~ ko) 



A A 2 \2{n-k )+l 

at oo is normalized. This completes the first step of expanding the series (12.11) into a continued 
fraction (see [24]). 

As was shown in the set of the normal indices of coincides with the following 
sequence 

n 2 — k < ■ ■ ■ < nj — k < . . . . 

Now one can apply the above reasoning to the function ipx and so on. By recursion we obtain 
the following P-fraction 



e eoEtbl £j~i £ j b2 j 



(2.4) 



3»j-l 



Po(A) - Pi(X) - ■■■ ~ Pj(X) 



where Ej = ±1, bj > and pj(X) = X kj +p^LxX kj ~ 1 + ■ • -+p^ X+p ^ are real monic polynomials 
of degree kj (see also [TT] . [14]). Note, that rij = k + fci + • • • + kj-x- 

It also should be mentioned that there exist functions for which the set {bj,p^\ . . . ,p^-i '■ 
j G Z + } of coefficients of the P-fraction is not necessarily bounded. In particular, the Cauchy 
transform of the signed measure constructed in [30] gives such an example. 
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The continued fraction (12. 4p can be considered as a sequence of the linear-fractional 
transformations (see [21], Section 5.2]) 



T» : - 



Pj(X) +ejtiju 
having the following matrix representation 

(2-5) W,(A) = ^ Y j e Z+. 

The superposition To o Ti o • • • o Tj of the linear-fractional transformations corresponds to the 
product of the matrices W;(A) 

(2.6) W [0J] (A) = K?(A))J fe=1 := W (A)W 1 (A) . . . W,(A). 

To give an explicit formula for W[o,j] in terms of pj, bj, Ej, define the polynomials P, + i(A), 
Qj+i(\) by the equalities 

The relation W [0 j](A) = W[ Q , i _i](A)W i (A) (see (|2^]l ) yields 

(2.8) W M (A) (J)= W[0J _ 1|( A) ( ^ ) = ( ) , i 6 N. 

So, the matrix W[o,j](A) has the form 

< 2 - 9 ' >w*> - ( 7v?w 1&$ ) • ^ z +- 

Further, the equality 

~p£<$ ) = >W W) ( ? ) = £>W*) ( 'fa ) ■ i € n, 

shows that the polynomials P/(A), Qj(X) are solutions of the difference equation 

(2.10) Ej-i£jbj-iUj-i — pj(X)Uj + bjUj +1 = (j G N), 
obeying the initial conditions 

Po(A) = l, P 1 (A) = ^, 

(2.11) 60 

Qo(A) = 0, Qi(A) = ^. 

According to (12.91) . the (J + l)-th convergent of the continued fraction (12.41) is equal to 

h := T o T x o ■ ■ ■ o T, (0) = -Q j+1 (X) / P j+1 (X) . 
The relations ( \2.9L ( 12. 6ft . and (12.51) imply the following statement. 

Proposition 2.1 Q14J). T/ie polynomials Pj, Qj satisfy the following generalized Liouville- 
Ostrogradsky formula 

(2.12) eMQi+i (A)Pi(A) - g i (A)P i+1 (A)) = 1 (jG Z + ). 
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3. Generalized Jacobi matrices 

The main goal of this section is to present a special class of generalized Jacobi matrices. 

Let p(A) = p n X n + ■ ■ -+pi\+po be a monic scalar real polynomial of degree n, i.e. p n = 1. 
Let us associate to the polynomial p its symmetrizator E p and let the companion matrix C p be 
given by 

fO ... ~Po\ 
1 o - Pl 





fPi ■ ■ 


■ Pn\ 




(3.1) 


E P =i\ . 


" oj 


C p — 




\Pn 







1 -Pn-lJ 



As is known, det(A — C p ) = p(X) and the spectrum c(C p ) of the companion matrix C p is simple. 
The matrices E p and C p are related by (see [20] ) 

(3.2) C P E P = E p C p . 

So, CpEp is a symmetric matrix. 

Definition 3.1 ([13], [23]). Let pj be real monic polynomials of degree kj 

Pj (X) = \ k > + pgL^" 1 + ■ • ■ + rf ] A + p£\ 
and Zet = ±1, 6j > (j G N). T/ie tridiagonal block matrix 

(A q Bo 0\ 

Bo A 1 B 1 

3.3 H = 

1 ; Bi A 2 ■■ 

V° '•/ 

where Aj = C p . and kj + i x fcj matrices Bj and kj x matrices Bj are given by 
(3-D Ojf= ( -1,^ 



6, 



.0 ... 



,o 







toiZZ 6e called a generalized Jacobi matrix associated with the sequences of polynomials {£jPj}^ 
and numbers {bj} c . 



I oo 
tj=0- 



Remark 3.2. The papers [13] . [m| . and [23J are only concerned with the case of generalized 
Jacobi matrices which are finite rank perturbations of classical Jacobi matrices. In fact, the 
generalized Jacobi matrix in question is associated to the P-fraction (12 .4p or, equivalently, to 
the sequence of matrices VVj having the form (12.51) . 

Prom now on, we suppose that 

(A.l) there exists N E N : degp.,- < N, j E Z + , 

(A. 2) swp{bj, \p^\ : j E Z+, k = 0, . . . , % - 1} < +oo. 



Let if Qoo ^ denote the Hilbert space of complex square summable sequences (w ,Wi, . . . ) 



with the usual inner product. Setting 



(3.5) 



n = 0, nj ^J2 k i 0' gN )' 



i=0 
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define a standard basis in £? ^ by the equalities 

e i,fc = {^,n j +Jfe}~ (J e Z + ; A; = 0, 1), e := e ,o- 
Define the symmetric matrix G by the equality 

(3.6) G = diag(G'o,G 1 ,...), C, (j G Z+). 

Further, we may identify via usual matrix product the matrix G with an operator on the 
linear space £o of finite sequences of £f 0oo y Its closure will be also denoted by G. In view 
of ( 1A.1I) . ( 1A.2I) . the operator G defined on £f 0oo \ is bounded and self-adjoint. Moreover, G~ l is 
a bounded linear operator in £^ n. 

Let Hy t q be a submatrix of if (G), corresponding to the basis vectors {%k\i=j' \ <_1 

(0 < j < / < +oo). The matrix iiyn will be called a finite generalized Jacobi matrix. 

Let 53[o,oo) be a space of elements of ^f 0oo ) provided with the following indefinite inner 
product 

(3-7) [x, y] = (Gx, y)e [Q aa) ( x , v e 4>oo))- 

Let us recall [3] that a pair (.£), [•, •]) consisting of a Hilbert space $) and a sesquilinear 
form [•, •] on Sj x is called a space with indefinite inner product. A space with indefinite 
metric (#,[•,•]) is called a Krein space if the indefinite scalar product [•, •] can be represented 
as follows 

[x,y] = (Jx,y)sj x,yefi, 
where the linear operator J satisfies the following conditions 

J = ,r l = J*. 

The operator J is called the fundamental symmetry. So, one can see that the space fj[o,oo) is 
the Krein space with the fundamental symmetry J = signG (see [3] for details). 

Proposition 3.3. Under the assumptions flA.ll) . flA.21) . the considered generalized Jacobi ma- 
trix defines a bounded self-adjoint operator H (a generalized Jacobi operator) in the Krein space 
53[o,oo), that is, 

(3.8) [Hx, y] = [x, Hy] x, y G % )0O )- 

Proof. It is not hard to see that, according to flA.ll) . (IA.2p . the matrix in question generates a 
bounded operator in £j[o )0 o)- Relation (13. 8p is implied by (13. 2p (see [13] for details). □ 

Let us extend the system {-P ) (A)}°f =0 , {Q.,'(A)}^L by the equalities 

(3.9) P hk (X) = X k P J (X), Qj,k{X) = X k Qj(X) (j eZ +] k = 0,...,k J -l). 
Setting 

P M (A) = (P li0 (X), JW-i(A), ■ ■ ■ . ^-o(A), . . . , P i)% -i(A)), 
Q[u](A) = (Qi,oW, Qim-iW, Qj,oW, Qj,k 3 -i(X)). 

one can rewrite the system (12.101) . (I3.9P in the following manner 

(3.10) P m (X)(X-H [0tj] ) = (0,...,0 J b J P J+lfi (X)) (j G Z + ), 

(3.11) Q m {X){X-H m ) = ( 0,...,0,- £ Q ,0,...,0,& j g j+1 ,o(A)) (j G Z+). 

Since Qo,o(A) = ■ ■ • = Qo,fc -i(A) = 0, the relation (13.111) reduces to 

(3.12) Q {ljj] (X)(X-H lhj] ) = (0, . . . , 0, bjQ j+ i j0 (X)) (j G N), 

It follows from (13.101) and (13 . 1 21) that the eigenvalues of if[o,j] and iip,j] coincide with the roots 
of Pj + \(X) and Q J+ i(A), respectively. 
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Proposition 3.4 ([13]). The polynomials Pj and Qj (j £ N) can be found by the formulas 

(3.13) Pj(X) = (6o...fei-i)- 1 det(A- J ff [0)i - 1 ]), 

(3.14) Q 3 (X) = e (6 ...^-i)" 1 det(A- Pf^-i]). 

The formulas (13. 13[) and (13.141) in the classical case can be found in (HI Section 7.1.2] and [H 
Section 6.1]. The following statement is an easy consequence of the recurrence relations ( 12.101) . 

Proposition 3.5 ([13].[14]). Let j G N. Then 

i) The polynomials Pj and Pj+i have no common zeros. 

ii) The polynomials Qj and Qj+i have no common zeros. 

iii) The polynomials Pj and Qj have no common zeros. 

Taking into account the equality Gioj]H[o,j] = HZ^Gpj] which is implied by (13.21) (see [13] ) 
and setting 

n [o,j]W = G [oj] P [o,j]( A ) T ' £[o,j](A) = G [oj]Q[o,j](A) T , 
one can rewrite (I3.10p . (13.111) in the form 

(3.15) (A - F [0|i ])7r[ Q)i ](A) = jh,!). :.a{X)< j {) (j G Z+), 

(3.16) (A - iy [OJ ])C[oj](A) + e , = e i 6 i Q i+ i,o(A)e J - >0 (j e Z+). 
Further, let us set 

(3.17) tt(A) = G- 1 (P ,o(A), . . . , P ,*o-i(A), • • • ) T , £(A) = G- 1 (g ,o(A), . . . , Qo,fco-i(A), • • • ) T - 
Now, it follows from (13 .151) - (13. 171) that the following formal equalities hold true 

(3.18) (A - H)n(X) = 0, (A - H)£(\) = -e , . 

The first equality in (13.181) allows us to characterize the point spectrum of H. 
Proposition 3.6. A G & P (H) if and only if tt(X) G ^? 0oo \- 

By the definition of 7r(A) and the assumptions (lA.ip . f]A.2[) . we see that 

oo 

vt(A)g4 oo) ^^|P :; -(A)| 2 <+oo. 

4. Weyl solutions and Weyl functions 

If for some A G C there exists a solution {Wj(X)}'jL Q of the recurrence relations (I2.10p 
such that 

(4-1) W(A)}£o G 4 >oo) and {Wj(X)}™ + {^(A)}^ 

then we will say that there exists a Weyl solution {W / J (A)}°2 =0 of the recurrence relations (I2.10p 
at the point A. Since {Pj(A)}°^ and {Qj(A)}°^ are linearly independent solutions of (12.101) . 
the Weyl solution admits the following representation 

(4.2) W J (X)=Q 1 (X)+m(X)P J (X), 

where m(A) is a complex number. The following statement shows the relation between m(A) 
and the operator H. 

Proposition 4.1. Let X G p(H) and let 

(4.3) m(X) = [(H-X)- 1 e,e], e := e 0fi . 

Then the family {Wj(X)}jZ given by (14.21) is the Weyl solution of (12.101) at X. Moreover, 
{Wj(X)}f =0 has the form 

(4.4) W j {X) = [{H-X)- 1 e,e j)0 \, j G Z + . 
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Proof. For A G p(H) the relation (13.181) implies that there exists a number m(A) G C such that 

e(A)+m(A) 7 r(A) = (iJ-A)- 1 eG^ oo) . 
Since G is a bounded operator, we have 

(4.5) G(£(A)+m(A)7r(A))G4 iQo) . 
Using (13.171) yields 

\{H - \y l e, e jfi ] = (G(£(A) + m(A)vr(A)), e ij0 )^ = Q 3 (X) + m(X)P j (X) =: W 3 {X) (^ P,(A)). 

The latter relation means that {Wj(X)} ( jL is a solution of (12.101) . Due to (14.51) . we obtain that 
{Wj{X)}f =Q G £f 0:Oo) . Now, it follows from flXTH that 

W (A) = m(A) = [(# - A) _1 e, e}. 

□ 

Definition 4.2. The function m defined by (14.21) is called a Weyl function of the operator H . 

Remark 4.3. It should be mentioned that a general treatment of the Weyl functions of clas- 
sical Jacobi matrices in the framework of extension theory of non-densely defined symmetric 
operators was proposed in [26J. A general treatment of the Weyl functions of symmetric oper- 
ators in Krein spaces in the framework of extension theory was presented in [T5]. The function 
m defined by ( 14. 3 p on p(H) is also called the m-function of H (see [13] . |18j). 

Since H is bounded, m admits the representation 

oo 

(4.6) m(X) = [(H-X)- 1 e,e] = -J2^ (\M > \M) 

3=0 

where = [H z e, e]. 

Analogously, one can define m-functions of shortened generalized Jacobi matrices. 
Definition 4.4. The function 

(4.7) m m {X) = [{H M -X)- l e,e] 
is called the m-function of i?roji ■ 

Making use of the structure of H[qj], we obtain that (see [T3] ) 

/ x / x det(A - #rwi) 

4-8 m m (X) = -e -Ml . 

det(A - H[oj]) 

According to (I3.13p . (13.141) . the formula (14.81) can be rewritten as follows 
(4-9) m m (X) = 

It follows from ( I2.10p (see [13] for details) that the m-function m[oj](X) and the m-function 
m[ij](A) of H[i t j] are related by the equality 

(4 ' 10) ^"' tfl^W 3eK 

An analogous statement for infinite generalized Jacobi matrices is an essential ingredient 
in the proof of the following result. 

Theorem 4.5. Under the assumptions ( lA.ip . ( 1A.2[) . the generalized Jacobi matrix H is uniquely 
determined by its Weyl function m. 
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Proof. By using the Frobenius formula, one can see that the Weyl function m of H and the 
Weyl function 

m[i j0O )(A) = [(#[i,oo) - A) _1 ei i0 ,ei i0 ] 
of #[i,oo) are related by the equality 

(4.H) m W = n m + 7vm — m' W > ll^ll > ll^ooll 

(a more detailed reasoning can be found in [13]). Further, consecutive applications of the 
relation (14. lip leads to the P- fraction ( 12.41) . So, one can uniquely recover the generalized 
Jacobi matrix H. □ 

Remark 4.6. In the definite case, formulas (14. 8ft and (14. 9ft are easy consequences of the theory 
developed in [26]. In this case, it is well known that the Weyl function determines the classical 
Jacobi matrix uniquely (for instance, see [TS], [2S])- It is worth to mention that this result 
as well as Borg type uniqueness result was recently extended to the case of normal matrices 
(see [27]). Besides, a canonical form of a normal matrix (an analog of Jacobi matrix for 
selfadjoint matrices) was also introduced there. Some inverse problems for finite generalized 
Jacobi matrices were considered in [12], IT3l. 



5. The resolvent set of H 

Here, following the scheme proposed in pQ, the characterization of resolvent sets p(H) 
of generalized Jacobi operators is obtained. We begin with an auxiliary lemma which gives a 
criterion of the density of ran (H — A) in £^ 

Lemma 5.1. For A G C the equation 

(H - \I)x = e j)k 

has a solution x = x(j, k) G £f n for all j, k G Z + iff there exists a Weyl solution at A. 

Proof. 1) First, let us consider the case where j — k — 0. It follows from (I3.18P that the 
equation (H — X)x = e ,o has a solution belonging to £\ ^ if and only if there exists a number 
m(A) G C such that £(A) + m(A)7r(A) G £f 0oo y Furthermore, in this case we have 

(5.1) x = x(0,0) =£(A) +m(A)7r(A). 

2) Next, let k = and j G N. According to (13TT5j) and (j3TTB|) . we see that 

(5.2) (H - A)£[ j](A) = e , - e j b j Q j+1}0 (X)e ji0 + e j b j Q jfi (X)e j+ i fi , 

(5.3) (H - A)7r [0)i ] (A) = -£ j 6 j P J+ i i0 (A)e J - o + ^I^JX)* s . Ul . 

Adding (I5.2p multiplied by — Pj and (15.31) multiplied by Qj, one obtains 
(5.4) 

(H-X) [-P,-(A)e [0j] (A) +Q j (X)7r m (X)] = -P i (A)e 0i0 -e i 6 J -(P i+1 (A)Q^(A)-Q i+1 (A)P i (A))e ij o. 
Due to (12.121) . (15 .4p can be rewritten as follows 

(H - A) [-P,-(A)^ [0>i] (A) + QiWnodV] + p jW^o = e if0 . 
The latter relation shows that the equation (H — \)x = ej t0 has a solution x = x(j, 0) G £? ^ 
if and only if eo,o £ ran (i7 — A) . Moreover, the solution admits the representation 

(5.5) x = x{j, 0) = -P,(A)£[ j](A) + Qj(X)iri 0t j](X) + P,(A) [^(A) + m(A)7r(A)] . 

3) Finally, assume that k ^ 0. Observe that 

(5.6) (P~ - A)e ji0 = e jtl - Xe jfi , ...,(H- A)e ijfe 2 = e ijfe x - Xe jjk 2 - 
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The chain of equalities (I5.6P imply that the equation (H—X)x = ej^ has a solution x = x(j, k) E 
^fo oo) ^ anc ^ on ly ^ e i,o ran (-^ ~~ A) ( or ) equivalently, eo,o G ran (if — A)). Besides, the solution 
x = x(j, k) can be expressed in the following manner 

x = x(j, k) = ej . fc _ 1+ Ae ijfc _ 2 +. • •+A fe (-P :; (A)e [ o, j] (A)+g,(A)7r [ o ii] (A)+P :; (A) [£(A) + m(A)vr(A)]). 

□ 

Remark 5.2. In fact, Lemma 15.11 gives a way to express the formal inverse operator Pa to 
H — A. Moreover, it is not so hard to see that the Weyl solution at A exists if and only if 
e G ran (H — A). 

Now we are ready to prove the main result of the present paper. 

Theorem 5.3. Under the assumptions flA.lj) . flA.21) . A G p{H) if and only if there exist a Weyl 
solution {Wj(\)}'jL Q at A and numbers q G (0, 1), C > such that 

(5.7) \PiWW j {\)\<Cq n i- ni , i<j. 

Proof. Let us prove the sufficiency Let TCk '■= spanje.,^1 ? G Z + } and let Vn k be the orthogonal 
projector onto Hf. in ^? 0oo v We start with proving the boundedness of the operator P A P Wfc for 
any k E {0, . . . , N — 1}. First, it is convenient to consider the operator GR\Pu - Taking into 
account (15.51) . we obtain 



(GR\Vn Gjfl)i,k — (GR\ejfl)i,k '■— (GPac^o, e;^) 



X k P i {X)(Q j {X)+m(X)P j (X)), i < j; 
X k P j (X)(Q i {X)+m(X)P i (X)), i>j. 



It is clear that one can represent the operator GR\Vh as the sum of upper and lower triangular 
operators: GR x V Ho = P A 1} + P A 2) . To be more precise, we choose R^ in the following way 

R\ e j,k = for A; ^ 0, P^ e J)0 = yf\ where 

(1) _jx k P i (X)(Q j (X)+m(X)P j (X)), i < j; 
[Vj kk ~\0, i>j. 

Now, one can prove that P A is bounded. Indeed, setting 

l kj-l l kj-l 

j=0 k=0 j=0 k=0 

we get the following relation 

l / l ki-l 

(R^x,y) e 2 = ^2x jfi I P^e^o,^^^,^ 

j=0 \ i=0 k=0 

I j ki — 1 

= E^'.o E E A fe P i (A)(Q,(A) + m{X)P 3 {X))y hk . 

j=0 i=0 k=0 

Thus, (15.71) yields 

i jki—i 

\(R^x,y) P \ < Cmax{l, lA^ 1 } ]T |x i>0 | £ £ q n ^ ni \y i>k \- 

j=0 i=0 k=0 
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Notice that rij — rii = k\ > ^ 1 = j — i and, therefore, we have 

l=i+l l=i+l 

l j h-i 

\(R^x,y) P \ <Cwi*{lA\\ N - 1 }Y t ^^^\xM\\Vi*\ 

j=0 i=Q k=0 
I l kjs-l 

< Cmax{l, {X^^j^q 8 ^ \x jfi \ \Vj-s,k\ 

s=Q j=s k=0 

< CViV max{l, lA^" 1 } Vg 4 ||x|| f2 ||y||, 2 = C 1 ~ ^ \\x\\ P \\y\\ P . 

Hence, R^ is a bounded operator. Similarly, one can prove the boundedness of R^ ■ So, we 
have proved that GR\P-h is bounded. Since G~ l is bounded, R\Vt-c is also bounded. 
Further, let 6 {1, . . . , iV — 1}. From (I5.6P one can deduce 

R\Zj,h — e j,k-l + ^R\Gj,k-l = VkGj,k + ^R\Vk^j,k, 

where V k : ej t k i— > &j t k-i is an isometric operator from Tik to 'Hk-i- If h G Tik then R\h = Vkh + 
\R\Vkh. So, the boundedness of V\ and R\Vn implies that R\Pm is bounded. Analogously, 

N-l 

RxPjik is bounded for k G {2, . . . , N — 1}. This implies that R\ = R\Pni is a bounded 

operator. The latter means that the domain of R\ is £f 0QO y Since A is not an eigenvalue, we 
have ker (if — A) = {0} and ran (H — A) = £? s. Now, applying the Banach theorem on inverse 
operators we obtain that A G p(H). 

The necessity of (15. 7p follows from [10] and the relation 



(R\ej,o)i,ki 



Pi(A)(Q i (A)+m(A)P i (A)), i < j; 
P J (A)(Q,(A)+m(A)P l (A)), i>j, 

which directly follows from (|5.5|) . □ 

Remark 5.4. In the case of nonsymmetric tridagonal operators Theorem 15.31 was proved in [1]. 
In fact, we have extended the scheme proposed in [1] to the case of generalized Jacobi matrices. 
A similar result for banded matrices with nonvanishing extreme diagonals in terms of the 
corresponding vector polynomials was obtained in [8]. 

6. The Floquet theory 

In the present section, by using Theorem 15.31 we give a description of spectra of periodic 
generalized Jacobi operators. 

Definition 6.1. Let s G N. A generalized Jacobi matrix satisfying the properties 
Ajs+k = ^k-, Bj s+ k = Bk, Sjs+k = £fc> j £ ^+? k G {0, . . . , s — 1} 

will be called an s-periodic generalized Jacobi matrix. The corresponding generalized Jacobi 
operator in 7i[o,oo) will be also called an s-periodic generalized Jacobi operator. 

Evidently, any s-periodic generalized Jacobi matrix satisfies the assumptions flA.ll) . flA.2j) 
and we have 

(6.1) W js+k (\)=W k (\), JGZ + , fcG{0,...,s-l}. 

The main tool for analysis of periodic generalized Jacobi operators is the following matrix 

- W[o,,-i](A) - I £s _ lbs _ lPs _ liX) Ps(X) 
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The matrix T(A) is called the monodromy matrix. Using (16. II) . we get the following relation 

(6.2) W [ o, JS+fc -i ] (A)=T^(A)W [ o, fe -i ] (A), j G Z+, fc G {1, . . . , s}. 

Let iui = wx(X) and w 2 = W2(A) be the roots of the characteristic equation det(T(A) — w) = 0. 
Introduce the following notations 

E := {A G C : \wi(X)\ = \w 2 (X)\}, E p := {A G C : P a _x(A) = 0, |6 s _ 1 g s _ 1 (A)| > |P(A)|}. 

Now, we are ready to give a description of spectra of periodic generalized Jacobi operators. 

Theorem 6.2. The spectrum of an s-periodic generalized Jacobi operator has the form 

a(H) = E U E p , a p {H) = E p . 

Proof. Since detT(A) = 1, we have that 

w 1 (X)w 2 (X) = 1. 

Step 1. First, let us prove that 

{A G C : P S _ X (A) ^ 0, \ Wl {X)\ ± \w 2 {X)\} C p(H). 
To be definite, assume that |u>i(A)| > |u> 2 (A)|. In this case, we see that 

(6.3) T(A) = h f) ( W f> w ° .) ( * 7 2 ) , det £ ^) = 1. 
Now, (I6.2p can be rewritten in the form 

Further, ( 16. 4ft is reduced to the following relation 

(a r\ iAi /\\ ( XxX^Wi — X2X3W2 —XiX 2 w{ + XiX 2 W J 2 \ /,n 

(D-oJ KV[ j S+ fc-i](AJ = j j j j yV[o,k-l]{X). 

\X3X4W1 — X4X3W2 —XsX 2 W 1 + X1X4W2J 

Multiplying (16.51) by the vector (1 0) T we obtain 

Sk-ih-iQja+h-iW \ = fx!X 4 w{ - x 2 x 3 w J 2 -Xix 2 w{ + XiX 2 w J 2 \ ( -e k - X b k _ x Q k _ x {X) 
£k-ibk-iPjs+k-i{X) J \x 3 x A w{ — x^x 3 w{ — x 3 x 2 w{ + xiX4w{) \ e fe _i6 fe _iP fe _i(A) 

Thus, for the polynomials Q. one has 

(6.6) Q js+ k-i{X) = -w^XxXiQu^X) + x^Ph-ii})) + w 2 (x2X3Qk-i(X) + xix 2 Pk-i(X))- 
Similarly, for the polynomials P. we have 

(6.7) P js +k-i(X) = -wKxzx^Qk-iiX) + x 3 x 2 Pfe-i(A)) + ^(x 4 x 3 Q fc _i(A) + x^P^iiX)). 

Note that x 3 7^ 0. Indeed, if £3 = then according to (16. 3p we would have that P s _i(A) = 0. 
Now, formulas (16. 6p and (16.71) yield 

(6.8) Q JS+ k-i(X) - -P js+k -i = -wi f-P fc _i(A) + Q k -i{X)) . 

x 3 \ x 3 J 

For brevity, define Cj.(X) := — ( — Pfc_i(A) + Qk-i(X) ) . Then the equality (16. 8ft can be rewrit- 
es J 

ten in the following way 

(6.9) Q js+k -i{X) - —Pjs+k-i = (w^-'C^X), w 2 = wl'% k G {1, . . . , s}. 

x 3 
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Since |iw 2 | < 1, it follows from (16. 9p that £(A) 7r(A) G ^f 0jOO )- A linear independence of 7r(A) 

and £(A), and (16.71) . (16.61) imply 7r(A) ^ ^foooV Now we are rea dy to verify the condition (15. 7p 
of Theorem 15.31 Let us assume that i < j. Then 



g,(A) - ^p,(A) 



|(-^A (fel) (A)+^/r ; (A))^C fe2 (A)|, ^ = V, fc^fe G {!,..., s}, 



where k\ = z(mods) and k 2 = j(mods). Since (rij — n^/N - 
j — i we have that (n^ — n^/N < j — i. Thus, one obtains 



E h)JN< 

l=i+l 



£ N \/N 

l=i+l J 



X3 



< |w 2 | j "^(A) < q n ^ ni C{X), 



where C(A) = sup {(|jf l} (A)| + |/ 2 lK2j (A)|)|C fc2 (A)|} and g = w 2 < 1. 

fc,fc 2 e{i,...,s} 
S^ep 2. Let us show that 

{A G C : P s _x(A) = 0, |6 S _ 1 Q S _ 1 (A)| > |P S (A)|} C a p (H). 

In this case the monodromy matrix can be represented as follows 



r(*a), 



T(A) 



-e a _i& a _iQ a _i(A) -Q a (A) 
P(A) 



1 Xi 

1 



wi(A) 
w 2 (A) 



1 



where u>i(A) = — e s _i6 s _iQ s _i(A) and u> 2 (A) = P S (A). Further, we have that 



w{ w{x\ — W 3 2 X\ 



-£fc-l&fe-l<5fe-i(A) 
£fc-i&fc-iPfc-i(A) 



-£fe-ifefc-iQjs+fe-i(A) 

£fc-l&Jfc-l-P/s+fc-l(A) y V W 2 

that is, the following formulas hold true 

Pjs+k-iW = w 2 Pfe-i(A), 
<5j S +fc-i(A) = -w{Q k -i(X) + (u^ - w^)xiP fc _ 1 (A). 
Since |u>2(A)| < |u>i(A)|, Gn(\) G $ 0oQ \ and, therefore, 7r(A) G ^f 0oo y So, we have proved that 



E p C <x(Pf). 

Step 3. By the same reasoning as in the Step 1 it can be shown that 

{A G C : P s _i(A) = 0, |6,_iQ 5 _i(A)| > |P S (A)|} c p{H). 

Step 4 ■ We complete this proof by proving the following inclusion 

{XeC:\w 1 (X)\ = \w 2 (X)\}Ca(H). 

Since Wiiu 2 = 1 we have that |u>i(A)| = |w 2 (A)| = 1. If T is diagonalizable then, due to (16.61) . 
(HTTP , there exist numbers (A) , (A) , 0jp (A) , (5^ (A) G C such that 

(6.10) Q JS+k - 1 (X) = w 2 J a < )} ) +wiaf, j G Z+, fc G {1, . . . , s}. 

(6.11) P is+fc _ 1 (A) = V'4 1) +^/3f, j'gZ + , fce{l,...,s}. 

Since the sequences Q. and P. are linearly independent, the sequence {Qj +mPj}™ is nonzero 
for any m G C. Also, observe that the sequence {lu^'a + (3\} ( jL l (a 7^ 0, \w\ = 1) converges only 
for w — 1. Taking into account these observations one can conclude that the series 

00 00 s — 1 

£ |Qj(A) + mP,(A)| 2 = E l«£° + + ^(4 2) + } )! 2 



j=0 



i=0 k=0 



does not converge for any number m G C. So, we have that £(A) + m7r(A) ^ 
m G C. Similarly, one can conclude that 7r(A) ^ ^f 0oo )- 



'[0,00) 



for any 
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If T is similar to a Jordan block then w\ = w 2 and Wiw 2 = 1. So, W\ = ±1. First, let us 
consider the case w± = 1. The monodromy matrix takes the form 



T(A) 

Further, (I6.2p yields 



Xi x 2 \ A 1 \ / x 4 -x 2 

x 3 X4I \0 1/ I— X 3 X\ 



yV [0js+k ^(X) = ( X2Xi ~ XlX f. ~ X2X * ~ Xl + ^ ~ XlX2 ) Wr^-xKA). 

Thus for any m G C the vectors ^(A) + m7r(A) and 7r(A) do not belong to (-f 0oo y Analogously, 
one can consider the case Wi — — 1. □ 



Remark 6.3. A description of spectra of classical Jacobi operators was obtained in [17] . 
Let us remind that Wi and w 2 are the roots of the equation 

(6.12) w 2 - (P S (X) - £ s ^b s Q s ^(X))w + 1 = 0. 
Due to (16.121) . one has 

(6.13) wi + w 2 = P S (X) - £ s _i6 s Q s _i(A), 

(6.14) Wl w 2 = 1. 

Remark 6.4 ([7]). Formulas (16.131) . (16. 14ft allow us to give another description of E 

(6.15) E = {X G C : (P S (A) - e 8 _ 1 6 a Q._ 1 (A)) G [-2, 2]}. 
Example 6.5. Let us consider the following difference equation 

(6.16) ^u j+1 - X\ - = (j G N). 

Clearly, the three-term recurrence relations H6. 161) generate the following 1-periodic generalized 
Jacobi matrix 

(A B \ 
E-[b A :.}. »). J" 

In this case, by easy calculations, we have 

P (A) = 1, Pi(A) = 2A 2 , Q (A) = 0, Qx(A) = 2. 

Since Po(X) = 1 7^ 0, it follows from Theorem 16.21 that cr p (H) = 0. Next, according to 
Theorem 16.21 and (16. 15f) , we have 

a{H) = {X G C : 2A 2 G [-2, 2]} = [-1, 1] U [-i, i]. 

By the same reasoning as in [7], one can prove the following statement on the structure 
of the set E. 

Proposition 6.6 (|7J). The compact set E has no interior points. The open set D := C\E is 

connected. The functions w% and w 2 are single-valued in D. 
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7. PADE APPROXIMANTS 

Our goal in this section is to prove convergence results for Pade approximants. 

+00 

Definition 7.1 Q28J). The [L/M] Pade approximant to the function (f(\) = — is 

j=0 

defined as a ratio 

Al L / M ] (I) 



BMW (i) 

of two polynomials A^ L I M \ B^ L ' M ^ of formal degree L and M , respectively, such that B^ L ^ M ^(Q) ^ 
and 

+00 

£ + / [L/M1 ( A ) = 0(A~ (i+M+1) ) (A - 00). 

In the case L = M = n, the [n/n] Pade approximant is also called the n-th diagonal Pade 
approximant. 

Let us consider the Weyl function m(X) = [(H — A) _1 e, e], where the corresponding matrix 
H satisfies the assumptions (1A.1I) . (1A.2I) . The representation (14.71) of mpj-i] yields 

(7.1) m|0 ,_ 1| (A) = -^ = - T ---...-^ rT + O^j (A — > 00), 

where Sj = [ff l e, e] (z = 0, . . . , 2rij — 2). Moreover, it was shown in [13] that m[ j_i](A) has the 
following asymptotic expansion 

2nj— 2+fcj / 1 \ 

(7-2) m [o,i-i] (A) = - J2 lfe + *S^*7 



8=0 



where Sj = [if l e, e] (i = 0, . . . , 2rij — 2 + fcj). The latter means that the rational function 



(7.3) /l"^l(A) = m [0 , ; _ 1 |(A) = A B[nj/nii ^ = f _ _ (J = 1,2,...) 



A^ A > 

is the [n_j/nj] Pade approximant to m. Besides, it follows from the Pade theorem (see [5j 
Theorem 1.4.3]), that for L and M satisfying 

L > Uj, M >nj, L + M < 2n,j + kj-l 

the [L/M] Pade approximants coincide with /KA^l, and for L, M satisfying 

L <nj + kj - 1, M < rij + kj - 1, L + M > 2uj + kj 

the [L/M] Pade approximants do not exist (for details see [13J). 
In what follows we need the following definition. 

Definition 7.2. The set 0(H) := {(Hy,y)p. : = 1} C C is called a numerical range of the 
operator H . 

Clearly, the numerical range of a bounded operator is a bounded set. By the Hausdorff 
theorem we have that a(H) C 6 (if) (see [22J). 

Theorem 7.3. Let m(X) = [(H — A) -1 e, e] be the Weyl function of H satisfying (lA.lj) . (1A.2j) . 

JTien i/iere exists a subsequence of diagonal Pade approximants f^ n ^ n ^ to m, which converges 
to m locally uniformly in C \ 0(H). 
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Proof. First, note that 0(fT[ O)T ,]) C ®(H), and, therefore, 0(if[o, n ]) C Q(H). As a consequence, 
we have that if A G C \ Q{H) then A G p{H %n] ) for all n G Z + . Let A G C \ Q{H). Then 
ran (H — A) = £? % and for any finite vector cf) we have 

(# [0>i] -A)-V->(#-A)-V, j — > +00. 

Thus, one can obtain 

(7 4) /KM](A) = mpy-nCA) = [(%,-!] - A) _1 e, e] 

= ((# [0J -i] - A) -1 e, Ge), 2 -> ((i? - A) -1 e, Ge)*. = m(A). 
According to [221 p. 336 Theorem 3.2] one has 

1 1 



dist(A,0(^ [Ol n])) " dist(A,0(#))' 



Hence, the family /K/ n jl = m^j] is uniformly bounded on compact sets in C\Q(H). Actually, 
the following estimates hold true 

II 1 1 

i[oj](A)| = \((H[o,j] ~ Xy 1 e,Ge) £ 2\ < \\(H [0J] - A) _1 e||^ \\Ge\\p < 



mil 



dist(A,0(#)) 

So, the family _/"[%'/%] i s uniformly bounded and, therefore, f^/ n ^ is precompact. Now, to 
complete the proof it is sufficient to apply (I7.4p and the Vitali theorem. □ 

Remark 7.4. In the proof of Theorem 17.31 we used the method proposed in [7J for complex 
Jacobi matrices. Note that Theorem 17.31 is a generalization of [T9| Theorem 1] . More precisely, 
we do not suppose existence of all diagonal Pade approximants and all poles of the existed 
diagonal Pade approximants belong to the convex set 0(H). 

Further, following the scheme proposed in [TJ, we find out a behavior of the associated 
polynomials at the points of the resolvent set. 

Proposition 7.5. For any A G p(H) the following inequality holds true 

(7.5) limsuplP^A)! 1 ^ > 1. 

Proof. In fact, the proof is in line with that in [TJ . However, we give the proof here for the 
convenience of the readers. Using (12.121) and (14.21) . one can see 

(7.6) P,(A)(5,6,H/ J+1 (A)) - P j+1 (A)( £ ^-(A)) = 1. 
Due to (IA~TD . f[A2D . and (15771) . one has 

(7.7) MiWi(A)| <C iqx \ 0< 91 <1. 

It follows from (17. 6p and (17. 7p that the sequence P/(A) can not be majorized by a geometric 
sequence pP with p < 1/qi, that is, for any p < \jq\ and any positive constant C 2 the inequality 

|Pi(A)| <<V 

is not satisfied for an infinite number of indices j. We thus have 

\P jk (\)\>C 2 pi*. 

If we choose p such that l/qi > p > 1 we obtain the required result. □ 

From these statements one can deduce the following result on convergence of diagonal 
Pade approximants. 

Theorem 7.6. Under the assumptions (lA.lj) . ([A. 21) . for any A G p(H) there exists a subse- 
quence of diagonal Pade approximants to m(A) = [(H — A) -1 e, e] ; which converges to m(X) at 
A. 
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Proof. From (17.51) we get that there exists a subsequence jk such that 
(7.8) \P h {\)\>C 2 p>\ P>1- 

Further, Theorem 15.31 and (I7.8P imply the relation 

W jk (X) 



q - 
P 



which completes the proof. □ 

Acknowledgments. I express my gratitude to Professor V.A. Derkach for reading the manu- 
script and giving helpful comments. I also thank the anonymous referee for helpful remarks. 
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